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Abstract 

We consider controlled differential equations and give new estimates 
for higher order Euler schemes. Our proofs are inspired by recent work 
of A. M. Davie who considers first and second order schemes. In order to 
implement the general case we make systematic use of geodesic approxi- 
mations in the free nilpotent group. 

As application, we can control moments of solutions to rough path 
differential equations (RDEs) driven by random rough paths with suffi- 
cient integrability and have a criteria for L 9 -convergence in the Universal 
Limit Theorem. We also obtain Azencott type estimates and asymptotic 
expansions for random RDE solution. When specialized to RDEs driven 
by Enhanced Brownian motion, we (mildly) improve classic estimates for 
diffusions in the small time limit. 

Introduction 

We consider controlled differential equations of the form 

dy = V x (y) dx 1 + ... + V d (y) dx d 

on the time interval [0, 1]. When the Revalued driving signal x and the vector 
fields are Lipschitz continuous then there exists a unique solution for every 
starting point yo- Building on recent work by A. M. Davie |7j, we control the 
Holder norm of y in terms of suited Holder norms of x and its iterated integrals 
and also obtain higher order Euler estimates. To be able to go beyond the first 
and second order case discussed by Davie we use ideas from sub-Riemannian 
geometry. En passant, we construct geodesic approximations which respect the 
geometry of the vector fields. Finally, by passing to the limit, all estimates 
extend to solutions of rough path differential equations (RDEs) in the sense of 
T. Lyons. 

As application, we show that random RDE solutions driven by a sufficiently 
integrable geometric rough path are in L q (fi) for all q < oo. The examples we 
have in mind arc RDEs driven by Enhanced Brownian motion (in which case 
we are effectively dealing with a Stratonovich stochastic differential equation) 
as well as certain Enhanced Gaussian - and Enhanced Markov processes. 
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When specialized to RDEs driven by Enhanced Brownian motion, the above 
mentioned Euler estimates are closely related to classic estimates by Azencott. 
In fact, at least in absence of a drift vector field, our main Euler estimate 
sharpens a key result in jj when applied to diffusions en temps petit, which is 
precisely the result which is quoted and used in Ben Arous' and, later, Castell's 
work HE]. 

As the reader may suspect, the robustness of the rough path approach allows 
to obtain Azencott type estimates for arbitrary random RDE solutions driven 
by a sufficiently integrable geometric rough path. The case of RDEs driven by 
EBM, as discussed in the last paragraph, is just a special case of this general 
result. 

Notation 1 The dimensions o/R d ,R e are fixed and will not appear explicitly 
when we write out the dependence of constants. In general, constants that appear 
in lemmas, propositions, theorems etc have an index that matches the number of 
the statement. In the proofs we indicate changing constants by a running upper 
index. 

1 Preliminaries: Euler scheme of order N 

Define = ©f =0 (R d )® k , by convention (R rf )®° = R. Let x be an Revalued 
Lipschitz path and define the k th iterated integrals of the path segment x|r Sj t] 
as 

ft fU k f'U2 

g**.-.<* := / / ... / <!.,-;, ..A,-'!' . 

J s J s J s 

and so that g fc = (g fc '* 1 '"'' l *) i i e ^ d j £ (R d )® . For later convenience 

set g° = 1 G (R d ) 8 = R. We then define the (step- AT) signature of the path 
segment x|r Sit ] as 

x^ftvCz^l + fy eTW (R d ). 
fc=i 

We say that a vector field is in Lip 7 (R e ) if it has |_7j bounded derivatives 
and the [7] * -derivative is {7}-H61der continuous. 

Definition 2 Given vector fields V\, V& & Lip 1 (R e ) (= bounded & Lipschitz 
continuous vector fields) and a M. d -valued Lipschitz path x on [0, 1] , we let y = 
7r (0,yo; x ) denote the unique solution to the (control) ODE 

d 

dy t =J2 Vi (»*) dx * = v fa) dxt > 1 e I°> ^ ■ 
1=1 

started at j/q- 
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The following lemma is left as a simple exercise. 

Lemma 3 Assume that {Vi) 1<i<d G Lip 1 (K e ) . Let x be an R d -valued Lipschitz 
path on [0, 1] and let y t — 7r (0, tjq; x) t . Then, for all < s < t < 1, 



where depends on (the Lipschitz norm of) the vector fields Vy, ...Vd- 

Let us now define the Euler approximation of order N to a control ODE of 
the above type. To this end, let H denote the identity function on R e and recall 
the identification of vector fields with first order differential operators. 

Definition 4 Given {V l ) 1<l<d G Lip^ (R e ) , g G T {N) (R d ) and y G K e we call 

N 

/" v «: / y,N,g] : E E 1, • • • 1,. ' ' " . 

e{i,. ..,<*} 

ifte (increment of) the step-N Euler scheme. 

This definition is explained by 

Lemma 5 Assume that (Vi) 1<i<d G Lip^ (R e ) . Lei x 6e an M. d -valued Lipschitz 
path on [0, 1] and let y t = it (0, yo', x) t . Then, for all < s < t < 1, 



E / • • • H (y ri )-V h ---V iN H (y s )} dx\\ ■ ■ ■ dx% 

■ __ J s<r 1 <...<r N <t 



tl ,...,in 

e{i,...,<0 



and there exists a constant depending on N and V\, Vd such that 

\y s>t - p.,N,x»*\ <Cfa(J \dx r 

Proof. Let / be smooth and note that Vi G Lip^ implies Vi 1 ■ ■ -V. k f is C 1 for 
1 < k < N. Iterated use of the fundamental theorem of calculus gives 

N-l 

f{Vt) = f(Vs)+J2 E / V tl ---V lk f(y s )dxl\---dxH 
fe=l i u -,ik Js <n< - < r k<t 
£{l,...,d} 

+ E / V il ---V iN f{y rx )dx^---dx%. 
e{i,. ..,<*} 
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This first part is then proved by specializing to / = H. For the second state- 
ment, lemma |21 gives 



\Vs,t 



< 



Gj^J^ \dx r 



Lip^-regularity of the vector fields implies that Vi 1 ..Vi N H (■) is Lipschitz 
and hence, for r £ [s,t], 

\V n ..V lN H(y r )-V n ..V lN H(y s )\<q^J^ \dx r \. 



This leads to 



t \ N +! 



[ [V ll ..V lN H{y ri )-V ll ..V lN H{y s )]dxl\---dxll <<%( I \dx r \) 

J s<n<...<r N <t \J s / 



and summing over the indices finishes the proof. 



2 Preliminaries II: Algebra of Iterated Integrals 

The set Tf (R d ) = {g G T N (R d ) : g° = l} is a group under truncated tensor 
multiplication: if g = 1 + g 1 + ... + g N = 1+ g and similar for h then for 
k = 0,...,N 

k 

\fc \ ^ „i -nx — i 



(g®h) fe = Y l g i ®b k 



i=0 



The neutral element ise = 1 = 1 + + .. . + and the inverse is given by the 
usual power series calculus 

(1+gr^i-g+r 52 -... 

For every A e R, the dilatation map S\ is defined componenwise by g fc i— > A fc g fc , 
k = 0,...,N. 

Sx : (g fe ) » (A fc g fc ) , A G K. 
Obviously, Tf (R d ) is a Lie group. Its Lie algebra can be identified with 

T N (R d ) ee {geT N (R d ) :g° = 0}, [g, h] =g®h-h®g 

and the exponential map with exp : Tf (R d ) -> Tf (R d ), g ^1 + g+ig® 2 + ... 
We recall some well-known facts. See CHI d] for further references. 

Proposition 6 (Chen, [15 ) Let x : [0, 1] — > R d be Lipschitz continuous with 
(step-N ) signatures x s . t = Sn (x) s t . Then 

S N {x). t ® Sn {xY = S N (x) s _ u . (1) 
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We define G N (R d ) = exp (L N (R d )) where 

L = L N (M d ) = R d © [R d ,R d ] © [R d , [R d ,R d ]] © ... c (R d ) , 

qN tod\ j s a Lj e subgroup of (R d ) with respect to ©-multiplication and 
known as step-N nilpotent free group over 18L d . 

Theorem 7 (Chow, |18|) For every g G G N (R d ) i/iere exists an R d -valued 
Lipschitz path x such that Sn (x) x = g. More precisely, G is the group 
generated by {exp («):t£ R d } so f/iaf every g G G is the signature of a (finite 
number of) concatenation of straight path segments. 

Theorem 8 (Geodesic Existence, |T5]) .For every g G G w (R d ), 

|g|| := inf | y |7 t | rft : 7 : [0, 1] — > R d Lipschitz continuous, 7 (0) = 0, Sat (7)0 1 = 8 
is finite and achieved at some minimizing Lipschitz continuous path 7*, i.e. 

||g|| = / and S N (7*) 1 = g. 

Jo 

Moreover, by simple reparametrization, we can state that for every s,t G R wf/i 
s < i i/ien exists a Lipschitz path x s,t : [s,i] — > R d wii/i signature g and length 

Wei- 

^(a;"'*) =g and / |cfcr s '* | = ||g|] . 

Remark 9 G w (R d ) can fee given a subriemannian structure so that the path 
t G [0,1] 1 — * S , at(7*) 04 is a subriemannian geodesic which connects the unit e 
with g G G N (R d ) , see 1 j<ffi . Thus, strictly speaking, 7* is not a geodesies but 
the projection of a geodesic. 

The geodesic existence theorem has useful consequences. Ifg, h G G N (R d ) then 
(i) ||g|| = iff g = e, (ii) symmetry: ||g|| = ||g -1 ||, (hi) sub-additivity 
||g®h|| < ||g|| + ||h|| and (iv) homogenity ||* A g|| = |A| ||g|| for all A G R, 
hold true. In particular, d (g, h) := ||g _1 © h|| defines a left-invariant metric on 
qN {^ d ^ t ne Carnot-Caratheodory metric. 

Theorem 10 ([18 ) (a) The topology induced by Carnot-Caratheodory metric 
coincides with the manifold topology of G N (R d ) and the trace topology as as 
subset of Tf (R d ) . 

(b) The map g i-> ||g|| is continuous in this topology. 

(c) The space G N (R d ) with metric d is Polish. 

Proposition 11 ([11 ) Let (i = 1,2) be continuous homogenous norms 

on G N (R d ) , that is, norms that satisfies properties (i) and (iv) and such that 
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g i— » |||g|||j is continuous w.r.t. r. T/ien i/iere exists a constant c E [l,oo) suc/i 
t/ia£ IH-lllj ~ |||-||| 2 by which we mean 

1 in in < in in < in in 

For instance, 

i til/* 

g = max g 

fe=l,...,jV' 1 

provides a useful example of a continuous homogenous norm on G N (R ) other 
than II -11. 



3 Preliminaries III: Geometric (Holder) Rough 
Paths 

Here, and in the remainder of this paper, we work exclusively with Holder 
modulus 

u (s, t) =t— s. 

Let p E [1, oo). A path x from [0, 1] to G N (R d ) is l/p-H61der continuous if for 
all s,te [0,1] 

||x M || <Cu;( S: t) 1/p 

for some constant C. This class is denoted by C' 1 / p - H61der ([0, 1] , G N (R d )). We 
can restrict attention to paths with pinned starting point. The (homogenous) 
l/p-H61der "norm" (there is no linear space here) on cVz-Holder ([0, 1] , G N (R d )) 
is defined by 

ll X lll/p-H61dcr:[0,l] = I l X H 1/p- Holder = SU P ' ~T7^ 

0<s<t<l LU (S,t) /P 

and there is a 1/p- Holder metric based on the CC-metric, 

«l/p-H61dcr:[0,l] l x , X) = "l/p-Holdcr ( x , X) = SUp rrj— . 

0<s<t<l Uj(s,t) IP 

We also set 

doo;[0,l] ( x > x ) = d oo (X, x) = SUp d (X 8) i, X s , t ) . 

0<s<t<l 

Theorem 12 ([10 ) (%) C 1 / p - mlder ([0, 1] , G w (R d )) is a complete metric space 
under the metric dx/ p -Holder- 

(it) Every l/p-Hdlder continuous path x E C 1 '^" 616 -^ ([0, 1] , G N (R d )) can be 
approximated by Lipschitz paths x n : [0, 1] — > R d m £/ie sense that 

Sn (xn) x uniformly on [0, 1] 



G 



and sup n ||5jv ( x n)\\i/ p .HSlder < 00 ■ ^ n / ac ^ we ccm /"^ Lipschitz paths x n such 
that 

Sup ||SW (sBn)|li /p .H fiWer < 3 IIxHj/ ffsWep . 
n 

(Hi) Assume p > 1. De/me C°^ p - HSlder ([0, 1] , G w (R d )) as the closure of lifted 
Lipschitz paths S N (x) under the metric d 1/p . H6lder . Forx £ C^ p - Hdlder ([0, 1] , G N (M d )) 
we have 

xeC °^P- mldEr ([0,l},G N (R d )) iffr(6;x)= sup l|XS|t J -> as J -> 0. 

0<s<t<l w (s,t) /p 

t — S<<5 

In particular, C°,VP-HMder Q 0) X ] ) G w ( R d)) c C l ' p - HSlder ([0, 1] , G^ (R d )) . 

One can see that C*°' 1 /P- H51dor ([0, 1] , G w (R d )) is Polish whereas (7 1 A'- HSlder ([0, 1] , G 
lacks separability. Recall that [p] denotes the integer part of some (positive) real 
number p. 

Definition 13 (pHHIi]) A path in C 1 / p ' mlder ([0, 1] , (R d )) is called a 
weak geometric p-rough path (with H older- control u> ) . A path in (jo,i/p-Hoider ^jq^ -q ^ q[p] 
is called a geometric p-rough path (with Holder- control u>). 

Proposition 14 (HS1Q7]) Let x € C 1 / p ' H5lder ([0, 1] , G^ (R rf )) and N > 
\p].Then there exists a unique lift of x to a G N (R d ) -valued 1/p-Holder contin- 
uous path w.r.t. Carnot-Caratheodory metric on G N (R d ), denoted by Sn (x) : 
[0, 1] — > G N (R d ). Moreover, there exists a const G (p,N) such that 

\\$N (x)\\l/ p -Hdlder;[nS] ^ ffifl I l X l I l/p-Hdlder;[0,l] • 

4 Generalized Davie Estimates 

In this section we show that the step-iV Euler approximation is a good approxi- 
mation to ODE solutions in small time, even if we control only the homogenous 
1/p- Holder norm of Sn (x) . In the case of N = 1,2 this result is due to A. M. 
Davie, [7]. The existence of geodesies associated to the Carnot-Caratheodory 
metric is our main tool to generalize his results to the step- TV case. 

Recall that a control ODE driven by Lip w vector fields has the step- TV Euler 
approximation 

7T (s, y s ; x) t » I y s ,N,S N (x) 6 



The Geodesic Existence theorem, applied to g = Sn (x) s t , yields the short- 
est path in M. d whose iterated integrals mimick the first N iterated integrals of 
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the path segement ac] • We called this path x s,t = (a;J'*)„ e r s t i • By construc- 
tion, its step- TV Euler approximation over [s, t] is exactly equal to I y <" N < SN ( x )s,t 
and we are led to the equally good step-N geodesic approximation 

ir(s,y s ;x) st « tt (s, y s ; x s ' t ) s t . 

This step- A'' approximation is sometimes easier to handle. It also respect the 
geometry given by the vector fields. Below, we shall use both. As last prepa- 
ration for the main result of this section, we need to understand the regularity 

y *-+l[y,N, g]. 

Lemma 15 Assume that (V % ) 1<l<d G Lip^ (R e ). For an element g G G N (R d ) , 

- 1*'"'*] < qm\y - y\ (ll«ll + Hell") 

where Cfjfil depends on N and the Lip N norm of the vector fields. 
Proof. By definition of the Euler approximation iy> N >& 

N 

iy ,N, s _jy,N, s = J2 £ [V il --V ih H{y)-V il ---V ih H{y)]^-' i K 

k=l ii,— ,ih 

e{i,-,<9 



Since y i— > • • -Vi k H (y) is Lipschitz, 



AT 



|^ s -^ g |<qfelE E If -SI I* 

fc=l ii,"' 

e{x,-,d} 

From equivalence of homogenous norms, |g fc,n '"' ,lk | < ^ ffg| llsll an( i hence 
l^.iV.g _ /fi ,Ar, g | <c^\ y -y\ max{||g|r , ||g||} . 

■ 

The next lemma is technical but very important. It quantifies the quality of 
step- AT Euler and geodesic approximations and gives ODE bounds which do not 
blow up with the Lipschitz norm of the driving signal. Recall that to (s,t) = t — s 
although the proof can adapted to general super additive control function |17j . 

Lemma 16 (Generalized Davie Lemma) Letp > 1 and {Vi) 1<i<d G Lip^ (R e ) 
for some integer N > p — 1. Assume that 

(i) x : [0, 1] — > M. d is a Lipschitz path with step-N lift x = Sn (x) and \\^-\\if p . Holder — 
Mi. 

(ii) yo £ R e with \yo\ < M2] Then, there exists a positive constant Cfjfil — 
G(Yfil (Mi), also dependent on p, N, M2 and the vector fields Vi,...,Vd but not 
dependent on the Lipschitz norm of x, such that for all < s < t < 1 

TT(0,y ;x) st <q Z g(M 1 )o; (2) 
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and with 9 = > 1, 



v 



Tr(0,y ;x) Sit -P°> N >s\<q m (M 1 )u J ( S ,t) 9 . (3) 
Moreover, if we assume 

(Hi) x s,t : [s,t] — » R d are Lipschitz paths such that Sn (x s ' t ) s t = x s ,t and such 
that 

J |dx s <*| < KM 1 uj(s,t) 1/p (4) 
for some positive real K then 

tt (0, y ; x) st - tt (a, y s ;x s ^) st < q^(Mi) w (a, t) e , (5) 



Cf7^| = Q^ti (Mi) a/so dependent on p,N,M2,K and the vector fields 



Remark 17 By hypothesis, ||x Si t|| < Milo (s,t) 1 ^ p and we can find a (pro- 
jected) geodesic 7*'[s, t] — > HH d with signature x Sjt and length J \dy*\ = ||x Sjt || < 

M\u) (s, t) 1 ^. In other words, paths {x s ^} as postulated in (Hi) always exist, 
even for K = 1 . ( By not fixing K we get a more general approximation result 
and highlight what is needed in the proof.) 

Proof. Without loss of generality, we assume that Mi > 1 (otherwise set 
Mi = 1.). 

Write y t = ir(0,y ;x) t and I\ f = y Syt - tt (s,y s ;x s ' t ) st . We first show © 
and divide the argument in two steps. 

First Step: Fix 0<s<t<u< 1. We try to control T s<u in terms of T s . t 
and T t .u- To this end, it is useful to define x s,t ' u to be the concatenation of x s,t 
and x*'". Observe that x s ' t,u \y s u i has the step-TV signature x Sit ® x t]U = x S) „ 
and 

\dx s ' t ' u \=J Idx^^+J \dx t ^ u \<2KM 1 uj{s,u) 1/p . (6) 
By uniqueness of ODE solutions, 

M ' Vs ' X >r-\7r (t, tt (s, y s ;x s <% ; x*>«) r if r e [t, u] . 

We have 

-r,,„ + r, lt + r t) „ = tt («, y 5 ; x s ' u ) Si „ - tt («, Vs - x s ^) s t - n (t, yt - x l > u ) t u 

= 7r(s,y s ;x s ' u ) su - ir (s,y s ;x s ' t > u ) su 

+ir (s, V.ix''*'") - tt (s, y s ; z 8 '*) - tt (i, y t ; ar*'*) 



By defintion of ir (s,y s ;x s,t ' u ), 

tt (s, y s ; a: 3 '*'") - tt (s, y s ; x"'*) = tt (t, tt (s, y s ; ar s '*) t ; a*'") 
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hence, 

-r s ,„ + r s , f +r t , u = n(s,y s ;x s ' u ) su -Tr(s,y s ;x s - t - u ) s u 

+n (i, tt (a, y s ; x s ' t ) t ; x^) tu ~ n (t, Vux^)^ . 
In particular, from lemma we have 

7r(s,y s ;x s '") su -7r(s,y s ;x s '*'") 

<\ir(s,y s ;x s > u ) S!U -iy*' N ^ + n (a, y s ; aj' Au ),, tt - 
< 2<flg(Mi) w (s, u)^ = 2^(^/0 u (a, u) 6 . 
using (gj and Then, 

< (tt (f, 7r (a, y s ; a^) t ; x^)^ - - 
, j7r(s,y,;a Slt ) t ,iV,x t; „ _ jj/ t ,JV,x tl . 



Once again, by lemma 

tt (t, 7r ( S , y s ; 7f\ ; s'")^ - f (^"%^ | + |^ ( t , y t ; x^)^ - J**- 

_ AT + 1 _ 

< 2Cgg|(M 1 )a;(i J «)— < 2Cjjg|(M 1 ) w (*,«)' 

Finally, by lemma IT51 

j7r(s,3/ a ;x s '*) t ,.!V",x t , u _ jy t ,N,x t 



< 



qjg|r,, t | max (m x w (t, , Mfw (t, 

< Cggir^l Mf max ( w (t, , u (t, u) N/p ) 

< 2Q^{M 1 )\T s , t \us{t,u) 1/p 
using that u)(t,u) < 1. Putting the pieces together, we have 

|-r s ,„ + r M + r t)U | < (2qJg(Mi) + 2qjg(Mi)) w (s, w) e +2qf S] (M 1 ) |r M | w (t, . 

It follows that 

|r,, u | < |-r, jU + r ait + r tiU | + |r, jt | + |r t , u | 

< \r a . t | (l + qJgtMi) « (i, U ) 1/p ) + |I\„| + qfgMr) w (a, M ) e (7) 
Second Step: For 0<s<i<u<l inequality J7J| can be rewritten as 



lu(s,u) u>(a,u) u> (s,t) 

uj(t,u) e \T tu \ 



uj(s,u) uj(t,u) 



10 



Define for r £ (0,1], 



g{r) 



sup 

0<s<t<l UJ (s, t)° 



Note that p (r) < oo. Indeed, this follows from 



|r. 



< 



y s ,t -Tr[s,y s ;x 

jV+1 



\t-s\ 



+ qj E pf 1 w(« l t) i!? 



The problem with this bound is that it blows up with The argument 

which follows shows that, in fact, p (r) will not blow up with |a;| iip ■ Pick arbi- 
trary points s < u such that w (s, u) < r, amd set t = (s + u) /2 so that 

u(s,i) = cj(*,u) = -w(s,u). 
We obtain from inequality J7J that 

lr -" < (D'em^i + cfaiMty* 



u> (s, u) 



< 



=7(1-) 



and taking the supremum over all s < u with u) (s, u) < r gives 

?M<2 1 -^(r/2) 7 (r)+qf S j 
After n iterated uses of the inequality for p we find 

n 

* • 
+qf i (M 1 ) 



fe=0 



fc-1 



fc=0 

A: 



3=0 



Let Cjgfr, Mr) := f^ 1 "^)* rfc 1 7 (&)) • Note that rR=o 7 (r/2 & ) 
is increasing in n and since 7 (r) < pu.G& Ml ^ r , the supremum over r of the 
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infinite product Jl^Lo 7 ( r 72 fe ) ^ s nn ite, which implies that sup 0<r<1 £ffffl( r ; Mi) 
is also finite. 



Hence, 



Q{r)<G^{M){2'- e ) n+l e / ' 



2 «+iy 1 1-2 1 - 6 



and sending n — » oo leaves us with (note > 1 here) 

e(r) < 

From the very definition of p with r = 1 we obtain 

ir^qJgCMiMM)' 



1-2 1 " 9 



TTwrd Step: Using lemma |21 and (@J, 



|dx s '*| 



< 



Then, for all s, t £ [0, 1] 



|2/ a ,t| < Vs,t ~ 7T (s,y a \x s>t ) s t + it (s,y s ;x s > t ) s t 

< C^(M 1 )Lj(s,t) e + C^CM l u>{s,tf /p 

< c^(M 1 )w(s,t) 1/p + c^M 1 w(s,t) 1/p 



= (fy(M 1 ) l j(s,t) 



Up 



Although C jj^| manifestly depends on X, we may specialize the construction 
using geodesies {x s, *}for which K = 1. With such paths, Gfei(Mi) of course 
would not depend on K. In particular, the Holder norm on y does not depend 
on K. 

Fourth Step: Finally, 10 is obtained from (J5J via triangle inequality and lemma 
taking into account (@J. ■ 

Corollary 18 There exists a constant Cff^ which may depend on p, N, Mi and 
the vector fields Vi,...,Vd so that for all Mi > 1 



This implies (the O -notation being understood as M\ — > oo) 

\nq m {M l ) = o({\nM l f). 
The same estimates holds for Gjfjfy allowing for additional dependence on K . 
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Proof. Inspection of the first step in the proof of Davie's lemma shows that 

qfa(^) . • > qfaic^o < qfei< +1 - 

The only difficulty is to control 

fe-i 



00/ j fc— 1 

< E (2^) n(i+qtHF" +l2 " 



fc=o \ i=o 

To understand the dependence of the right hand side on Mi, we define the 
function for some fixed a 6 (0, 1), 

A(fc,6) = a k U)zl (1 + &2-*') 
r(n,6) = ]TA(fc,6). 

We need to understand the dependence of hin^^oo T (n, b) on b. One could use 
a naive approach (the one used in the previous proof) to get 

A (fc, b) < a'TTto 1 (1 + &2-W) (8) 



TOO — 1 

^'=0 



< a^^ 1 exp (62- pj ') 



< a k exp 



1 - 2~f 



and hence, 

lim T(n,b) < — !— exp 

1 — a \ 1 — 2 p 

Unfortunately, the right hand side in the last equation grows to fast in b for our 
purposes. To obtain a better estimate, we first observe that 

A (fc, 0) = a k and - A (fc, b) = I £ j A (fc, 6) . 

Then we note that a: i-> 2" px / (1 + b2~P x ) = f 1 - (1 + b2-P x )~ 1 ) is de- 
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creasing in x so that 



k 1 2 _M 1 2~ p i 



^ 1 + 62-PJ ~ 1 + b ^ 1 + 

J=0 3=1 



< 



1 f°° 2~P x dx 



1 In (1 + 6) 
1 + 6 + 61n(2P) 



In 6 
* 3 X 



for all b > e. We also note that 



/ (b) := exp Q (In bf^j solves ^/ (6) = 3^/ (6) . 

and any other solution to this ODE must be a multiple of /. By ODE comparison 
we see that, for b > e, 

A(k,b) f(b) §(lnb)3 
A(fc,e)-/(e)" 

which implies that, using (jSJ) 

A(fc,6) < A(fc,e)e^ lnf ') 2 



JV+1\ 



< a fc exp(^ T -^ F j e l^^ 
After summing over all non-negative integers k we see that 

Um T(n,b)< C -^-^tll e i^)\ 

n—ioc 1 — a 

Hence, we have proved that 

exp ( ,_%_ v ) /o 
qSg(l,Mi) < _\ ^/ exp^-ln(q| H1 M 1 f 

< qf S] exp(3(^ + l) 2 ln(M 1 ) 

< qfg ] exp(6iV 2 ln(M 1 ) 2 ) . 

This lead to 

q^(Mi) < qf a Af 1 A ' +1 exp(6^ 2 ln(M 1 ) i 

< Cjfgpp (l27V 2 In (Aii) 2 ) , for M : > 3. 



By increasing Qfwi if needed we can assume that this estimate holds for all 
Mi > 1. Clearly, trie same estimate holds for ffifi|(Afi) • ■ 
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5 Euler Estimates for Rough Differential Equa- 
tions (RDEs) 

We consider controlled differential equations in the sense of T. Lyons. The 
driving signal is assumed to be a weak geometric p-rough path with Holder 
control w(s,t) = t — s. Recall that this means x : [0,1] -> GM (R d ) is 1/p- 
Holder continuous w.r.t. Carnot-Caratheodory metric on C?M (R d ). Lyons' 
theory |151 1171 ITC] then implies existence and uniquess of a solution to the 
differential equations driven by x along vector fields Vi,...Vd £ Lip p+e (R e ) 
started at some point yo £ R e at time 0. This RDE solution is also a (weak) 
geometric p-rough path, over R e instead of R d , denoted by 

7T (0,y 0) x) = y, 

with the same modulus of continuity as x. For our application it will be sufficient 
to consider the pathlevel RDE solution (obtained by projection) 

7r(0,|to,x) = y : [0,1] R e . 

Thus, y is anR e -valued 1/p- Holder continuous path in the usual sense. 

Theorem 19 Fix an integer N > [p] + 1 and Lip N -vector fields Vx, on 
R e . Let x be a weak geometric p-rough path with ||x|| u < M . Then there exists 
a unique pathlevel RDE solution 7r(0, j/0: x ) = V- Moreover, (a) there exists 
constant G[jg] = GfY^ (M), also dependent on p,N,yo and V\, Vd, such that 

\y\l/p-Hdlder;[0,l] — ( ~tlM 

and (b) a constant Grrn = QnT\{M) with similar dependencies such that for all 
0<s <t<l, 



y s ,Ar,5 JV (x) st ]| <q^{s,t) 9 with0=^j±>l. 



Vs,t - I 

Finally, keeping all parameters but M fixed, 

lnqTTjlnq^- O ((lnM) 2 ) as M -> oo. 

Proof. Lip^ +1 -regularity is more than enough to ensure existence and unique- 
ness of RDE solutions, see El GUI ■ 

(a) From Theorem 1 121 we can find Lipschitz paths x n such that 

S]p] (x n ) - x 

uniformly on [0, 1], such that 

ll 5 bl (^)H 1/p -H 5 lde r ^ 3M = M l 
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The Universal Limit Theorem implies a fortcriori that 

7r (0, y ; x n ) -> 7r(0,y ;x) 
uniformly on [0, 1]. On the other hand, Davie's lemma implies that 
sup|7r(0,y ;£ n )| 1/p _ H61dcr < q^i< oo 

n — 

where Cjj-g|is the constant Cj^|= Cfrfl|(Mi) from lemma [TBI It follows that 

k(o,y ;x)| 1/p _ H6kicr < q^|< oo. 

From corollary 1 181 

ln<jg|= O ((InAf!) 2 ) = O ((InAf) 2 ) . 

(b) By lemma ITU a weak geometric p-rough path x with Hx^ /p_Hoider — ^ 
lifts uniquely to a path S N (x)ex: [0, 1] -> G N (R d ) such that 



l X lll/p-H61dcr — 



for some constant Cj^g] — A). As in part (a) we can find Lipschitz paths 

x n such that 

S N (x n ) -> x 

uniformly on [0, 1], such that 

\\SN(x n )\\ 1/p . Hme r<^q^ = Ml 

Note that, by projection, (x n ) — > x uniformly on [0, 1] with uniform ho- 
mogenous 1/p- Holder bounds. As before, the Universal Limit Theorem implies 
that 

7r (0, yo; x 11 ) — > 7r (0, yo; x ) uniformly on [0, 1] 

while Davie's lemma implies the existence of Cj^Tft de facto Q^g] = ffig] (Mi) 
from lemma [TBI such that, uniformly over n, and for all < s < t < 1, 

tt (0, y ; x n ) s<t - I [tt (0, y ; z n ) s , A, ^ (jg»),J | < Cggp (*, t) fl 

with = (A + 1) /p. By continuity of the map 

0, g) G K e x G N (R d ) i > 7 [z, A, g] G M e , 
we can send n — > oo to obtain 

(0, y ; x) M - J [tt (0, y ; x) s , A, x Sii ] | < qfy (s, t) e . 
Finally, as above, 



In Cgg = O ((hiA/i) 2 ) = O ((In A/) 
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6 Asymptotic Expansions for RDE Flows 



We now consider RDEs driven by a random geometric p-rough path x = x (lu) 
defined on some complete probability space (f2, T, P) . We shall assume that the 



l/p-H61der;[0,l] 



has Gauss tails since this is the case for all examples we 



have in mind: Enhanced Brownian motion B, see [13 IE], Enhanced Fractional 
Brownian Motion B ff and other Enhanced Gaussian processes EI an d En- 
hanced Markov processes with uniformly elliptic generated in divergence form 
JT3J- However, the proof of the following theorem will make clear that the method 
works whenever the real-valued r.v. ||x|| 1 / p _ HS j t j er .r ji has some exponential tail 
decay. 

Theorem 20 Letp G (2,3), N > [p] + l and consider the random RDE solution 
7r(0,?/o! X ) driven by the random geometric p-rough path x (along fixed Lip N 
vector fields V\, ...,Vd), assuming that |] x |lx/p-ijoider'[o l] ^as Gauss tails, 



3a > : E 



expla \\X\\l/ p -Holder;[Q,l] 



< 00. 



The 



k(°>fo,x)| 1/p ._ ff5Her;[0)1] eL q (CL) for allqe [l,oo). 
Moreover, the remainder of the step-N Euler approximation is bounded in prob- 



ability. More precisely, there is a constant 
Vi, Vd such that for R > 1 and all t G (0, 1]^ 



dependent on a,p,N,yo and 



sup 

0<s<t 



7r (0,?/o,x) 



I 



y ,N,S N (x) J | > Rt^\ < q^exp (-e (haR)1/2/C M 



In particular, the l.h.s. tends to zero uniformly over t G (0, 1] as R —> oo anc 
the convergence is faster than any power of 1/R. 



Proof. By assumption, Mi = max |||x|| 1/p _ H61dCT . [01] , 
E [exp (aMf)] < oo. From Theorem IT§1 



l| has Gauss tails and 



sup 

0<s<t 



where 



tt (0, yo, x) 0)S - Iy , N , S N (x) 0;S | < q m {M l ) t^ +l ^. 

qtsi(Mi)<tioF (c fei lnMl)2 - 
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Therefore, 



sup 

0<s<t 



< 



7T(0,y Q ^) Qs ~P a ' N ' S ^o,. 
M > exp 



N+l 

> Rt~ 



Vin(fl/qy N 



< E I exp (aM 2 )] exp 



-a exp 



' 2y/]n(R/C^ \ 



< 



C^exp (_ e 0»(«)) 1/a /qft| 



where the last estimate is valid for every i? > 1 by choosing Ogn sufficiently 
large. ■ 



y ,N, S N (x) &?/7r (0, 7/0,2° 



Remark 21 27ie same result holds if we replace I 
where x 0,s is a geodesic associated to the element Sjy (x) s of G 

Remark 22 Even in the case of Enhanced Brownian motion, x = B, probabil- 
ity estimates of the unrestricted event 



sup 

0<s<t 



t(0,2/0)X) - I y ,N,S N (x) _ 



JV + 1 

> Rt~ 



valid for all (t, R) s (0, 1] x [1, oo) appear novel compared to the results given in 

0HOF 

The perhaps strongest estimate that has been extracted from Azencott's 
work in this context (see [U p 235]) is the following: in our notation (recall that 
7r(0, 2/o,B) solves a Stratonovich stochastic differential equation): 3a, c > : 
Vi? > : 



t->o 



limP sup 



0<s<t 



JV + l . 

> Rt~ < ce " 



(9) 



(Note that the exponent of t is(N + 1) /2 in contrast to (N + 1) jp in Theorem 
I20|) . We now show how © can be deduced from our general results. 

Proposition 23 We keep all assumptions of the preceding theorem but drive 
the RDE with Enhanced Brownian motion x = B. Then 0) holds with a = 
2/ (N + 1) and c — C fe| depending on N, y and V\, Vd- 



Proof. Choose p — p (N) s.t. 



iV + 2 
2< P <2—< 3 
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Then there exists a = a (p) > s.t. E [exp(aM 2 )] < oo where M = 

max {ll B lli/p-Hoidcr;[o,i] > 1 }- We set B = S n+i (B), well-defined by Proposi- 
tion [HI Then 



7r(0,y ,B) 0>s -P°^( B ^ 
From Theorem 1191 and Proposition 1 14| 



< 



p ,N+l,-B , s _ jy ,N,S N (x) 



7r(0,y ,x) 0)S -/ 3 ' ' JV+1 ' So -| < ^ A,+2 )/"e( c tisi lnM ) 2 
On the other hand 

jyo,N+l,S N +i(B) , B _ pa,N,S N (x) , e 



£ V il ---V iN+1 H(y )B* 



N+l,i lt — ,i N+ i 



n,...,*jv+i 
e{i,...,d} 



B 



W+l,tl,"- ,ijv + l 

0,s 



< 



Oral II Bo,. 



iJV+l 



Trivially, sup < s < t s^ 2 )/? = £( w + 2 )/p an d we are led to 
7r(0,2/ ,x) -/^< S »W°. 



sup 

se[o,t] 



< (l) + (2) 



iV + 1 



l ^ 2 + qf a su P ||B 0)S f +1 >^ 

s£[0,t] 



where with v = (N + 2) /p - (N + 1) /2 > 

(i) = p^q^p( c &i lnM ) 2 > 

and 



2 . i? 
2" 



Ar+l 



(2)=P| q|j U- 1/2 sup ||B , 



> 



R 



Gauss tails of M are more than enough to asset that (1) tends to zero as t — > 0. 
As for (2), Brownian scaling shows that (2) is in fact independent of t and hence 
equal to 



JV+l 



Cfei sup B 0iS 
U[o,i] 



*2 ' 



4 i/fiV+l 



> 



R 
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since Bup, 6[0il] ||Bo,.|| < ||B|| 1/p _ H61dCT . [01] < qijl|B|| 1/p _ H61der . [0)1] < q|jM\ 
Gauss tails of M now easly the claimed tail decay. ■ 

Remark 24 Estimate remains valid when I yo, N, Sjv (B) is replaced 
by 7r (0, yo] x°' s ) where x°' s is a geodesic associated to Sn (x) s . 

Remark 25 Note that a = 2/ (N + 1) is what one expects from integrability of 
the (N + 1) multiple Wiener-Ito integral which dominates the remainder as 
t — ► 0. Theorem \2(A gives an estimate valid uniformly int G (0, 1] and we would 
not expect an exponential tail. In fact, the given estimate implies a tail decay 
which is about as good as one can hope in absence of exponential decay and there 
seems little room for improvement. 

Remark 26 Only in absence of a drift vector field Vq does our step-N Euler 

approximation I yo, N, Sn (x) s coincide precisely with the approximation of 

£0 prop 4-3] and p 235]. When Vq ^ 0, the methodology and results of (p, q)- 
rough paths \14Jj could be used to produce the same approximations as those in 
the above cited references. If one accepts a few additional terms in the step-N 
approximation, it may be simplest to deal with Vq ^ via an RDE driven by the 
canonically defined time-space rough path. 

Remark 27 Proposition^^ is readily adapted to EFBM"B H with H > 1/4 and 
gives 



limP ( sup 

\0<s<t 



n(0,y ,B H ) 



yo,N,S N (H H ) o 



> Rt H ^ < ce 



Remark 28 In comparison to the Taylor expansion of Azencott, the approxi- 
mations by Ben Arous and Castell respect the geometry of the problem. This is 
also the case for our geodesic approximations although the efficient approxima- 
tions of sub-Riemannian geodesies remains a numerical challenge. 

7 L q Convergence in the Universal Limit Theo- 
rem 

We now give a criterion for L 9 -convergence in the Universal Limit Theorem. 

Proposition 29 Assume that for a random sequence of rough path 

swp\\x n \\ 1/p _ Hm [Q1] =M 

n 

has a tail Gaus tail. Assume that di/ p _jioider ( x n , x ) —* in probabililty. Then 
\n{0,y o ,x n ) -7T (0,yo, x)| 1/p . H6 . Wer . [01] -> m L q Vg e [l,oc). 
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Proof. By the universal limit theorem, 

Z n = |tt (0, ?/o, x„) - vr (0, yo,x)| 1/p . HSlder . [01] ^ in probability. 
At least along a subsequence, Z nk — * a.s. as k — > oo and 

k(0,y o ,x)| 1/p . HS i der;[O)1] = Um k(0,y ,x„J| 1/p . H61der . [01] 

and therefore 

Z„ < 2sup|7r(0,yo,x„)| 1/p _ H51dCT . [01] < 2C* 

n 

where C = C(M) G Vq 6 [l,oo) using Theorem I2U1 To obtain L« con- 
vergence we have to check that {Z^ : n > 1} is uniformly integrable. But this 
follows immediately from Z/-boundedness, r > 1, indeed 

E(|Z«| r ) < 2 9 ' r E(C(Af)« r ) < oo. 



Example 30 Let x = B be Enhanced Brownian motion over M. d with 

\Ml/p-Holder;[Q,l] < 00 where P G ( 2 ' 3 ) • 

Let (D n ) be a nested family of dissections of [0,1] with mesh \D n \ — ► and 
define the lifted piecewise linear approximation 

x n = S 2 (E [B\a [B t : t G £>„)]) . 

TTien observe that the area process of x n is obtained by the Levy area process 
condtioned on a (B t : £ € -D n ). Then 



sup||x„|| 1/p . K5Wer . [0)1] < supE 

n n 

Note that 

M (n) = E 
< E 



\l/p-H6lder;[0,l] 



a (B t : t G D n ) = sup M (n) = M < oo. 



Il/p-i?67der,-[0,l] 



U/p-HoWer;[0,l] 
and /rom £>oo6 's maximal inequality, 

q 



\a {Bf.te D n ) 
\o {Bf.tE D n ) 



1/9 



E 



Sup ||Xn|li/ p ./rfi(der;[0 1 l] 



< 



< 



1 



supE(M (n) q 



E 



ll/p-.f£o'«er;[0,l] 



Noting that (q/ (q — l)) 9 stays bounded as q — ► oo (in fact, it converges to e), we 
conclude that Gauss tails of Mi imply Gauss tails of M . Finally, fix p G (2,3) 
and repeat the above argument for p G (2,p) replacing p. By interpolation 

dl/p-Holderi^n,*-) ~ * a.S. 

and uniform Gauss tails follow from the trivial estimate || x n|li/ p -i2o7<ZeW0 1] — 

l X "lll/p-ffo7(ier-;[0,l]' 
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Remark 31 The same argument works for other Enhanced Gaussian processes 
with martingale approximations such as fractional Brownian motion, see ]1CH/ 
for the case H > 1/3. 
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